Partial Survival and Crossing Statistics for a Diffusing Particle in a Transverse Shear 

Flow 



o 
o 

(N 
D 

l> 
(N 



I 

o 
o 



> 
o 

o 
o 



X3 
O 

o 



X 



Alan J. Bray-'^ and Satya N. Majumdar^ 

^School of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, UK 
^Laboratoire de Physique Theorique et Modeles Statistiques (UMR 8626 du CNRS), 
Universite Paris-Sud, Bat. 100, 91405 Orsay Cedex, France 
(Dated: February 6, 2008) 

We consider a non-Gaussian stochastic process where a particle diffuses in the y-direction, dy/dt — 
rj(t), subject to a transverse shear flow in the x-direction, dx/dt — f{y). Absorption with probability 
p occurs at each crossing of the line x — 0. We treat the class of models defined by f{y) = ±v±(±y)" 
where the upper (lower) sign refers to y > {y < 0). We show that the particle survives up to time 
t with probability Q{t) ~ g^^j dgj-ive an explicit expression for 9{p) in terms of a and the 

ratio v+/v-. From 9(p) we deduce the mean and variance of the density of crossings of the line 
X = Q for this class of non-Gaussian processes. 

PACS numbers: 02.50.-r, 05.40.-a 



There has been a resurgence of interest in first-passage 
problems in recent years []|, particularly in the context 
of systems with many degrees of freedom Q . There are, 
however, relatively few exactly solved models involving 
non-Gaussian processes. This paper is devoted a class of 
such models, involving just two degrees of freedom, for 
which exact results can be obtained. These models are 
most simply described in terms of a particle moving in 
the two-dimensional plane {x,y), with stochastic motion 
in the y-direction and 'deterministic' motion in the x- 
direction (deterministic in the sense that the velocity in 
the x-direction depends only on the y-coordinate, which 
is itself however a stochastic variable). 

The class of models is defined by the following equa- 
tions: 
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where the upper (lower) sign refers toy > {y < 0), dots 
indicate time derivatives, and r](t) is Gaussian white noise 
with mean zero and correlator {ri{t)r]{t')) = 2Dd{t — t'). 
These models are non-Gaussian except for the case a — 1, 
V, which reduces to the random acceleration 
process x — V'q[t). 

Previous work on this class of models has addressed 
the first-passage problem in which the line x = is an 
absorbing boundary and the particle starts in the half- 
plane a; > 0. It was shown |^ that the probability, Q{t), 
that the particle survives until time t decays as Q{t) ~ 
for large t, with 
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where 
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The exponent 9 takes the value 1/4 for all a for the 'anti- 
symmetric' case 7 = 1 (for whichf(y) is an odd function) 
which has a simple explanation |j,|5j in terms the Sparre 
Anderson theorem . 

In the present paper we consider the case where the 
line a; = is a partially absorbing boundary, at which 
the particle is absorbed with probability p at each cross- 
ing. We have shown elsewhere that for such models 
the survival probability decays at late times as a power 
law with an exponent 9 that depends, in general, on p: 
Q{t) ^ t^^^PK Explicit analytical results for 9{p) have 
been obtained for a number of simple (and not-so-simple) 
models 7]. We also showed that the function 9{p) deter- 
mines the moments of the number of crossings of the 
boundary x = for the case where there is no absorp- 
tion. In the present work we use the methods of ref. j3] to 
determine the function 9{p) and the method outlined in 
ref. to investigate the crossing statistics for the prob- 
lem without absorption, obtaining closed form results for 
the mean and variance of the number of crossings in a 
given (large) time interval. These results for the mean 
and variance are, to our knowledge, the first such results 
for non-Gaussian processes. 

The derivation of 9{p) follows closely that of 9{0) in 
Ref. [3, so we only give the main steps. The principal 
difference from the former work lies in the boundary con- 
ditions imposed by the partially absorbing boundary. 

From Eqs. and Q we can write down the backward 
Fokker-Planck equation 

dQ 
dt 



(5) 



where Q{x,y,t) is the probability that the particle still 
survives at time t given that it started at position {x, y). 
The partially absorbing boundary at .t = implies the 
boundary conditions 

Q{Q+,-y,t) = pQ(0+,y,t), y > (6) 
Q(0+,-j/,i) = pQ{Q+,y,t), y > 0, (7) 
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where Q{x,y,t) is the survival probability for a model 
in which and V- are interchanged. It is clear that 
Q and Q are described by the same value of 9. In 
fact Q{x,y,t) = Q{—x,—y,t), since after interchang- 
ing v+ and V- the system is restored to its original 
configuration after a rotation by tt about an axis per- 
pendicular to the xy plane. The initial condition is 
Q(x,j/,0) = l = Q(x,y,0). 

Solving the full initial value problem is difficult so we 
will follow the approach used in ref. ^ by specializing 
to the late-time scaling regime where Q{x, y, t) ~ t^^'^P\ 
This approach exploits a generalization of the method 
introduced by Burkhardt [Sj for a — \ — ^ (the random 
acceleration problem). The idea is to extract explicitly 
the time-dependence expected at large t. This gives, 
asymptotically, 
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and a similar expression for Q{x, y, t), with scaling func- 
tions F±. In Eq. ijSJ, the functions F±{z) are the scaling 
function for y > and y < (— ) respectively. The 
(dimensional) prefactors (for y > and y < 0) in Eq. 
have been omitted since Eq. (01 is linear. The functions 
F+(z) and F-{z) are defined such that the prcfactor is 
the same for y > and y < 0. 

Inserting the form (jH)) into the backward Fokker-Planck 
equation (jSJ, we see immediately that the term dQ/dt 
leads to a term of order t^^^^^\ which is subdominant 
for large t and can therefore be dropped. The remaining 
terms give 

zF4'(z) + (l-/3-/32z)^^4(z) + 2/330F±(z) = O . (9) 

Expressed in terms of the variable u = P'^z, this equa- 
tion becomes Kummer's equation. Independent solutions 
are the confluent hypergeometric functions M{—2fi9, 1 — 
/3,/3^z) and U {-2139, 1 - 13,0^ z) The function 

M{—2(39, 1 — P,P'^z) diverges exponentially for z ^ oo, 
so must be rejected for y > 0. Thus we write 

W+/32yl/'^\ 



F+{z) ^ AU (^-2/30,1-/3, 
F_(z) = BU (-2139,1- (3, 
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with similar equations for F±{z) involving amplitudes A, 
B and C. 

Relations between the coefficients A, B, C, and the 
corresponding tilded variables can be obtained by impos- 
ing the boundary conditions lO and jT)), and requiring 
continuity of Q, dQ/dy, Q and dQ/dy at y = 0. After 



some straightforward but lengthy algebra, these bound- 
ary and continuity conditions eventually lead to a con- 
sistency condition on 9, which determines 9 — 9{p) as 
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where 



S = 



2p2cos2 +2sinh2 (iln- 

cos(7r/3) -I- cosh(ln7) 



(13) 



and we recall that (3, 7 are defined by Eq. Q). 

The general result, Eq. H12|l can be checked in a 
number of special cases: 

{i) p = 1. In this case S — 1 and 6* = 0. This is clearly 
correct since for p = 1 there is no absorbtion and 
Q{x,y,t) = 1. 

(ii) p = 0. For this case, 6 = (y^ - l/V7)^/[7 + 1/7 + 
2cos(7r/3)]. Inserting this into ifT^ . and carrying out 
some elementary manipulations, one recovers Eq. ©, as 
required. 

(iii) 7 = 1. This gives 5 = p^, and 
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The special case, (3 = 1/3, of this result corresponds to 
the random acceleration problem (a = 1) with partial 
absorption, and gives 9 = 1/4— (3/27r) sin~^(p/2), re- 
covering the results of Burkhardt |S| and De Smedt et 
aL^Qj for this special case. In the following, we exploit 
the general result H12I) to obtain results for the crossing 
statistics of the process defined by Eqs. ^ and 

In order to place the discussion of crossing statistics 
for non- Gaussian processes in the proper perspective, we 
begin with a short discussion of crossing statistics for 
Gaussian processes. Within the class of models discussed 
in this paper, only the case a = 1 = 7, corresponding to 
the random acceleration process, is Gaussian. 

For illustrative purposes, we begin with the random 
acceleration process, x — y, y — r/{t) or, equivalently, 
X = r]{t), with {ri{t)r]{t')) = 5{t — t'). For convenience 
we take the initial condition x(0) = = i;(0). Then 

x{t) = J* dt' /g dt" r](t"), and the two-time correlator is 
C{ti,t2) =< xlti){x{t2) >= i2ii/2-i?/6, where we have 
taken t2 > without loss of generality. The normalized 
two-time correlator is 



C{h,t2 
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for t2 > ti. Notice that C{ti,t2) depends on ti,t2 
only through the ratio ti/t2- This is an immediate 
consequence of the fact that, with the initial condition 
cc(0) = = a;(0), there is no timescale in the problem, 
so dimensionless correlation functions can only depend 
on ratios of times. If the initial position and velocity are 
non-zero, the scaling form (|15|l is recovered in the limit 
that ti and t2 are both taken large with the ratio held 
fixed. It should be noted that this scahng property of 
temporal correlations is not restricted to Gaussian pro- 
cesses, but holds for all the models discussed here. 

This scaling property implies that, if one introduces a 
logarithmic time variable [ill , T = In t, normalized 
correlation functions can only depend on differences of 
T-variables, i.e. the processes become stationary in log- 
arithmic time. Again, this applies to both the Gaussian 
and the non-Gaussian processes that we consider. For 
the random acceleration process, it follows immediately 
from (|15|l that the normalized correlator in logarithmic 
time, /(T) = {X{T)X[0)) where X{T) = x(t) / ^ {x^ {t) , 
is 

/(T) - (3/2) exp(-T/2) - (1/2) cxp(-3T/2) (16) 

for T > (and /(-T) = /(T)). 

For Gaussian processes there is a major simplification: 
the two-time correlator implicitly determines all proper- 
ties of the system, including the first-passage exponent 9 
and the crossing statistics. As an example, consider the 
mean density, p, of zero crossings in logarithmic time. 
The mean density p is finite only for smooth Gaussian 
processes whose correlator f{T) has the short time be- 
havior, f(T) = 1 - aT^ as T ^ where a = -/"(0)/2 
is finite. For such processes, the mean number of zero- 
crossings of the process X{T) in time interval T is 



dT'{5{X{r))\X{r)\) ^T{5iXiT'))) {\X{T')\) 

(17) 

since X{T) and X{T) are uncorrelated for smooth pro- 
cesses, {X{T)X{T)) = /'(O) = 0. Furthermore, X and 
X are normally distributed with variances /(O) = 1 and 
— /"(O) respectively. It follows that the mean crossing 
density is given by 



(18) 



a result first derived by Rice |l3l |. 

In a similar way it is possible to obtain an (albeit 
rather more complicated) expression for the variance, 
(n^) — {n)'^, of the number of crossings in time T, again 
expressed in terms of the correlator f{T), a result due to 
Bendat 0. 

For non-Gaussian processes, this approach no longer 
works. The first line, (n) = T{6{X) \X)\) stiU holds, but 
since the stationary distribution of (A", A) is in general 



not known, further progress seems to be impossible. We 
will show, however, that exact results can be obtained for 
the class of non-Gaussian models discussed in this paper 
by exploiting our knowledge of the function 9(p). The 
connection between the two was first noted in ref.0. 

Working in logarithmic time, as discussed above, we 
can write the survival probability Q(T) for the partial 
survival problem in the form 



g(r) = ^p"F„(T) 



(19) 



where P„ (T) is the probability of the process crossing the 
line X = n times in the (logarithmic) time interval T, 
and is the probability of surviving all n crossings. Now 
Q{T) decays asymptotically as Q{T) ^ exp(— 0(p)T), 
and the right-hand side of Eq. H19() can be written in 
terms of the cumulants of ri, to give 




exp(-0(p)r). 



(20) 



for large T, where {n^)c is the r'^ cumulant of n, and 
therefore 



f;(!£f(,0.^-«(p)T. 



(21) 



r=0 



for large T. 

The cumulants of n can now be determined by ex- 
panding both sides of Eq. (|21|l around p = 1, by writing 
p — 1 — e and equating coefficients of powers of e, to ob- 
tain the cumulants of n in terms of the derivatives of d{p) 
evaluated at p = 1. In this way one obtains 



[fi) = 
in) = 



-e'{\)T 
-e"{i)T 



(22) 
(23) 



etc. This approach gives, for the mean crossing density, 



P = 



sin(7r/3) 



T 



1 



27r/3 cos(7r/3) + cosh(ln7) 



(24) 



For the special case (3 — 1/3, 7 = 1 that corresponds 
to the (Gaussian) random acceleration process, we can 
check Eq. H24() against the general result H18|l that holds 
for any Gaussian stationary process. We find that both 
expressions reduce to p = •\/3/27r for this case. 

In a similar way we can calculate the second cumulant 
of the number of crossings by expanding to second order 
in 1 — p. The result for (ri^)c = {n"^) — {n)"^ is 



T 



[2Fsin(7r/3) - sin(27r/3)]. 



where 



F = l/[cos(7r/3) -I- cosh(ln7)]. 



(25) 



(26) 
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Again, this result can be checked for the random accel- 
eration process /3 = 1/3,7 = 1, which is Gaussian. For 
this case, Eq. gives {v?)c/T = 2/7r\/3. On the other 
hand, the second cumulant can be calculated exactly for 
any Gaussian process using Bendat's formula. This gives 
the result 2/7rA/3, as expected. 

It is noteworthy that for the class of models discussed 
here, arbitrary cumulants of the crossing number are 
readily obtained by simply computing the appropriate 
derivatives of 0{p), evaluated dXp — 1. Thus a knowledge 
of 9(p) provides a powerful tool. Even for Gaussian pro- 
cesses, a general expression for the moments of n above 
the second is not available for models where 6{p) is not 
known, which is essentially all but a few special models 

. . . 

It is also interesting to compute the mean time in- 
tervals (logarithmic scale) l± between crossings that the 
particle spends on the positive (negative) side of the X 
axis. For the special case 7 = 1 (where v+ = v-), it is 
clear that l-^- — I- — 1/p with p being the mean den- 
sity of crossings given by Eq. (|24(l . However, for arbi- 
trary 7 7^ 1, Z+ is, in general, different from Let 
p± denote the mean number of crossings of X = from 
the right (left) of the Y axis. Clearly p+ = p^ = p/2. 
To calculate the mean intervals l± in the general case 
(7 7^ 1) one can proceed as follows. We consider the nor- 
malized process X{t) in the logarithmic time scale r so 
that it is stationary. Let us first define a new variable, 
the 'occupation time', that measures the fraction of time 
spent by the process X{t) above (below) the X axis, 
L± = i Jg 9[±X{t)) dr. Taking the average, and using 
the stationarity, one gets {L±) = {9{zLX)). However, it 
is evident that (L±) = p±l±. Hence we get an expression 
for the mean intervals 



and f{x) = outside. Using this result, one gets the 



exact distribution of x{t). 



i± = —{e{±x)). 
p± 



(27) 



For the case where one has the symmetry X —X, 
such as the case 7=1, one gets, using {6{X)) = 1/2, 
the expected result Z+ = L = 1/p. However, for 7 7^ 1, 
the exact knowledge of p from Eq. 1)24(1 is not enough 
to calculate the mean size of intervals l± . One needs to 
compute, in addition, the quantity {9(±X)). To calculate 
this average we need to know the stationary probability 
density P(X), since {9{X)) = P{X)dX. 

The calculation of the probability density P{X), for 
7 7^ 1, is nontrivial. The only case where we have suc- 
ceeded in calculating P{X) exactly is the case where a = 
0. In this case, the equation of motion x{t) in the original 
time t reads, from Eq. i = v+9{y) — v-9(—y). Then 
one can write, x{t) — et + vTt where Tt — /q sign[y(t')] dt' 
is the sign-time of an ordinary Brownian motion, e = 
(v+—t;_)/2 measures the anisotropy and w = {v++v-)/2. 
The distribution of Tt for Brownian motion is well known 
to have the famous arcsine form of Levy Ts^, P{Tt, t) = 
where f{x) = l/[7r^/l~^] for^ -1 < a; < 1 



P{x.t) = JjG(^) where 
for — 1 -I- e/w < z < 1 + e/v 
Carrying out the integral of 



G{z) - l/[7T^l-iz-ey 
and G{z) = otherwise. 
P(x, t) only over the positive (negative) x axis, one gets 
{9{±x{t))) = ^± ^sin-^ (f). Using 9{X) = 9{x) and 
p = l/[7rcosh(ln7)] (obtained by putting /3 = 1/2 in Eq. 
(|24ll ) we finally get the exact mean intervals for the a = 
case: 



l± — 27rcosh(ln7) 



"1 

2± 






-sin-i 






TT 





(28) 



The determination of l± for other values of a (including 
the a — 1 case) remains an open problem. 

In this paper we have calculated the partial survival 
exponent, 9{p), for a class of (in general non-Gaussian) 
stochastic processes describing a random walker moving 
in a transverse "shear" flow. We have then used the re- 
sult for 9{p) to derive exact expressions for the first two 
cumulants of the crossing number (number of crossings of 
the line x = 0) working on a logarithmic timescale where 
the process is stationary. To our knowledge these are the 
first results for the statistics of crossing numbers for any 
non-Gaussian process. We have checked that our gen- 
eral result reduces, in special cases, to the known results 
for the random acceleration process. We have also com- 
puted exactly the mean time intervals between successive 
crossings of the y axis for the special case a — 0. The 
calculation of the mean time intervals for a 7^ remains 
a challenging open problem. 
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